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The  fracture  and  fatigue  studies  at  AMMRC  encompass  both  experimental  and  analy¬ 
tical  thrusts.  This  presentation  will  cover  a  few  representative  activities. 

The  status  of  test  program  on  fracture  behavior  of  unidirectionally  reinforced 
composites  using  double  cantilevered  type  specimens  will  be  reported.  Material  systems 
include  "S"  glass/epoxy  and  graphite/epoxy.  Limited  fracture  toughness  data  is  avail¬ 
able.  Additionally,  some  effort  has  been  initiated  on  cross  ply  specimens. 

Fatigue  studies  at  AMMRC  in  metals  reflect  a  continuing  interest  in  both  initia¬ 
tion  and  propagation  aspects.  A  technique  has  been  developed  for  detection  of  crack 
initiation  using  an  automated  photographic  process  for  recording,  periodically,  poten¬ 
tial  initiation  sites.  The  observations  are  discussed  in  conjunction  with  damage 
criteria.  Crack  propagation  studies  have  focussed  in  developing  a  "law"  which  accounts 
for  the  extremes  of  propagation  rates  at  threshold  levels  of  stress  intensity  and  at 
the  higher  levels  associated  with  unstable  growth.  Additionally,  crack  propagation 
experiments  are  described  which  are  aimed  at  exploring  the  transient  effects  on  crack 
propagation  of  sudden  changes  in  stress  intensity  levels. 

Supplementing  these  experimental  studies  are  extensive  analytical  efforts  directed 
toward  defining  the  stress  states  in  anisotropic  material  when  used  in  lap  md/or  mech¬ 
anical  joint  configurations.  Additionally,  extensive  effort  by  one  team  of  investiga¬ 
tors  has  been  devoted  to  analytical  techniques  for  crack  analysis.  These  ternniques 
have  been  applied  to  a  wide  variety  of  geometric  configurations  and  to  a  range  of 
material  types  including  anisotropic.  This  effort  will  be  described  briefly.  (Author) 
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I .  INTRODUCTION 


Within  the  Army  Materiel  Command,  AMMRC,  the  Army  Materials  and  Mech¬ 
anics  Research  Center  at  Watertown,  Mass.,  has  the  responsibility  for  research 
in  structural  materials  and  mechanics. 

The  scope  of  that  mechanics  responsibility  is  portrayed  in  Figure  1. 

I  shall  describe  the  effort  of  only  two  of  the  areas  shown  -  namely  fracture 
and  fatigue. 

In  the  fracture  regime  I  will  discuss  some  of  our  experimental  studies 
in  the  mechanical  behavior  of  layered  fibrous  composites  -  both  unidirec¬ 
tional  as  well  as  angle  ply;  I  will  mention  some  of  the  unique  characteris¬ 
tics  of  these  angle  ply  composites  in  producing  stress  concentrations  near 
free  edges  even  in  the  absence  of  gross  geometric  notches,  and  then  "'ll  cover 
our  purely  analytical  studies  in  fracture  mechanics,  particularly  as  oriented 
toward  anisotropic  materials.  Finally  I  will  conclude  the  fracture  dis¬ 
cussion  with  a  fleeting  reference  to  some  interests  in  dynamic  fracture. 

In  the  area  of  fatigue  I  will  describe  briefly  a  crack  initiation  study 
and  a  study  examining  the  transient  growth  rate  of  cracks  when  subjected  to 
a  step  change  in  the  applied  stress  intensity  factor. 

I  wish  to  irdicate  at  the  outset  that  my  role  as  speaker  is  as  a  re¬ 
porter  of  other  people's  work.  I  will  attempt  to  identify  and  acknowledge 
these  contributors  as  I  proceed. 
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II.  FRACTURE 


Slepetz  of  our  staff  is  carrying  out  an  in-depth  investigation  of 
the  fracture  behavior  of  unidirectional [1]  and  cross  ply  filamentary  layered 
composites.  Figures  2  and  3  show  superficially  the  nature  of  his  investiga¬ 
tion.  Using  specimens  of  the  compact  tension  type  shown,  he  has  obtained,  at 
least  for  unidirectionally  reinforced  composites,  load  deformation  curves  of 
the  type  shown.  These  characteristically  showed  a  series  of  successive  load 
drops  corresponding  to  "pop-ins"  early  in  the  loading  sequence,  but  these 
eventually  damped  out  and  disappeared  as  loading  and  crack  extension  proceeded. 
With  these  load  deformation  curves  and  experimentally  determined  compliance 
curves  as  shown,  the  end  result,  the  fracture  toughness,  was  easily  obtained 
using  the  relation 

Gc  =  (F2/2t)(dc/da)  (1) 

Representative  curves  for  two  material  systems  are  shown.  The  early  rise  in 
the  G£  vs  a  curves  is  associated  with  the  pop-in  behavior;  the  later  sharp 
rise  occurs  when  the  crack  approaches  the  far  end  of  the  specimen. 

Figure  4  shows  the  more  detailed  compliance  data  obtained  for  the  uni¬ 
directionally  reinforced  S-glass/ epoxy  specimen.  This  shows  compliance  C 
plotted  against  dimensionless  crack  length  a/w.  The  solid  line  represents 
the  equation  shown.  This  equation  is  based  upon  a  simple  strength  of  materials 
analysis  with  coefficient  determined  to  match  the  experimental  points.  The 
analytical  expression  was  then  used  to  determine  dc/da  for  use  in  Equation  1. 
Figure  5  shows  the  bulk  of  fracture  toughness  data  obtained  for  this  type 
specimen.  Note  as  in  Figure  2  the  early  pop-in  region,  followed  by  the  steady 
state  stage  and  the  ultimate  sharp  rise.  Figure  6  shows  the  corresponding 
toughness  data  for  unidirectional  reinforced  specimens  with  fibers  at  various 
other  orientations.  Again  the  same  trends  are  noted.  However,  it  is  noted 
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that  no  attempt  was  made  in  analyzing  these  latter  specimens  to  separate  ef¬ 
fects  of  the  crack  opening  stress  intensity  factor  from  the  in-plane  stress 
intensity  factor,  Kjj. 

The  next  figure  (Figure  7)  shows  some  of  the  effects  of  cross  ply  re¬ 
inforcement,  i.e.,  where  adjacent  layers  have  their  respective  fibers  in 
different  directions.  It  has  been  earlier  noted  that  for  unidirectional  com¬ 
posites  the  cracks  always  form  parallel  to  the  fibers;  in  cross  ply  cases 
the  cracks  still  tend  to  follow  the  fiber  direction  for  the  particular  layer 
in  which  the  crack  exists.  The  interaction  of  the  crack  systems  in  adjacent 
layers  leads  to  excessive  interlaminar  shear  stresses  and  ultimate  debonding 
or  delamination.  Furthermore,  as  the  "crack"  pi  gresses  by  this  mechanism, 
we  are  not  exactly  sure  of  the  appropriate  way  of  defining  a  fracture  tough¬ 
ness  quantity. 

To  get  a  better  base  line  for  understanding  the  behavior  of  these  fibrous 
composites  we  have  undertaken  an  examination  of  the  deformation  and  strain 
field  at  the  crack  tip  region.  Figure  8  shows  some  preliminary  results  of 
deformation  fields  in  a  unidirectional ly  reinforced  specimen  using  moire 
interference  techniques.  These  studies  were  carried  out  by  Sleprtz  and  Prof. 
Fu  Pen  Chiang  of  State  University  of  New  York  at  Stonybrook  College[2]. 

The  upper  row  of  interference  patterns  was  made  using  a  300  line/incli 
grating  and  white  diffuse  light.  The  second  and  third  rows  were  taken  using 
an  optical  multiplication  system  which  provided  600  and  1200  lines/inch  re¬ 
spectively.  The  columns  correspond  to  crossed  grids,  horizontal  grids  and 
vertical  grids,  respectively.  These  grid  arrangements  provide  all  the  defor¬ 
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mation  components.  Analysis  of  these  results  has  not  been  completed  yet. 

These  materials  characterization  studies  are  being  coordinated  with 
large  scale  panel  tests  incorporating  stiffeners  in  order  to  validate  the 


applicability  of  the  toughness  measurement  and  to  substantiate  recently  com¬ 
pleted  analysis  by  Lakshmikantham  which  will  be  discussed  later.  Specifically 
a  test  program  is  underway  to  observe  the  fracture  and  fracture  arrest  char¬ 
acteristics  of  stiffened  flat  panels.  Material  systems  being  examined  in  the 
panel  study  include  fiber  glass/epoxy,  graphite/epoxy,  and  boron/aluminum. 
Instrumentation  will  be  predominantly  strain  gage  and  moire  grids. 

Another  important  trait  of  angle  ply  layered  composites  is  the  edge  ef¬ 
fects  they  display;  significant  stress  concentrations  are  evidenced  even  in 
"simple"  unnotched  tension  strips.  Oplinger  et  al  of  our  staff  have  been 
investigating  this  both  experimentally  and  analytically[3] [4] .  Figure  9  shows 
the  trend  of  the  in-plane  deformation  field  as  well  as  the  interlaminar  shear 
stress  distribution.  Note  the  non-symmetry  of  the  deformations.  This  is 
caused  by  virtue  of  the  fact  that  the  filaments  in  each  layer  tend  to  re¬ 
orient  themselves,  under  load,  into  the  direction  of  the  tension  load.  In 
doing  this  there  is  a  scissoring  action  between  adjacent  crossed  layers  and 
this  sets  up  the  high  interlaminar  shear  stresses  shown.  Figure  10  shows  some 
experimental  results  again  using  moir£  gratings  to  determine  the  deformation 
fields.  The  upper  views  show  the  grating  direction  and  the  resulting  in-plane 
deformation  pattern  substantiating  the  predicted  trend  shown  in  Figure  9.  The 
lower  views  show  the  moir£  fringe  pattern  looking  in  on  the  edge.  The  quality 
of  the  fringe  pattern  is  not  good;  on  the  other  hand,  it  is  pointed  out  that 
this  specimen  is  only  .085"  thick,  and  techniques  for  applying  gratings 
under  such  circumstances  are  not  well  developed.  However,  again  the  fringe 
pattern  is  consistent  with  the  prediction  and  provides  substantiation  of  the 
presence  of  the  interlaminar  shear  stresses.  Figure  11  shows  the  variation 
of  the  fringe  pattern  with  applied  load  for  these  same  specimens. 
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The  effect  of  stacking  sequence  is  illustrated  in  Figure  12.  This  shows 
the  interlaminar  shear  distribution  in  a  24- layered  filamentary  composite  for 
two  stacking  sequences.  The  left  curve  is  for  a  clustered  sequence,  i.e., 
twelve  (12)  in  one  orientation  and  then  twelve  in  the  crossed  orientation. 

This  leads  to  a  rr  her  significant  stress  concentration  factor  of  approximately 
2.6  at  the  edge.  On  the  other  hand,  if  the  24  layers  are  stacked  in  alternating 
orientations,  i.e.  +6,  -0,  +0,  .....  then  the  peak  stresses  are  significantly 
reduced  and  the  stress  concentration  factor  is  approximately  1.27. 

Oplinger  has  also  treated [3]  the  problem  of  edge  effects  as  influenced  by 
the  structure  of  reinforcing  filaments  such  as  graphite  or  glass  yarns.  His 
results,  shown  schematically  in  Figure  15,  reveal,  for  example,  ;.hat  yarn  rein¬ 
forcement  gives  rise  to  more  pronounced  deformation  gradients  and  correspond¬ 
ingly  to  much  higher  interlaminar  shear  stresses  than  with  large  filament  re¬ 
inforcements  . 

Figure  14  summarizes  the  range  of  experiments  anticipated  in  connection 
with  these  edge  effect  studies.  Both  clustered  and  alternating  stacking  se¬ 
quences  will  be  investigated,  and  the  range  of  ply  angles  will  span  from 
+10°  to  ±45°.  Several  material  systems  will  be  explored,  i.e.,  boron/epoxy, 
graphite/epoxy,  and  boron/aluminum.  To  date  the  concentration  of  effort  has 
been  for  the  boron/epoxv  system  with  a  clustered  stacking  sequence  at  ±  20°. 

Oplinger  and  Gandhi  are  also  engaged  in  analytical  studies  of  joint 
design.  Two  generic  configurations  are  being  examined,  the  adhesively 
bonded  lap  joint  and  the  mechanical  (pinned)  joint.  As  suggested  in  Figure 
15,  the  focus  of  effort  here  is  on  layered  composite  constituents.  The 
layered  plate  with  circular  hole  is  being  examined  to  provide  input  for  the 
detailed  solution  of  the  pin  joint  problem.  As  a  preliminary  to  more  gen¬ 
eral  lap  joint  problems,  Gandhi [5]  has  completed  and  programmed  the  plane 
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strain  analysis  of  a  lap  joint  with  isotropic  constituents  and  zero  thickness 
bond.  The  configuration  and  results  for  both  the  bond  layer  shear  and  tensile 
stresses  are  shown  in  Figure  16.  The  singularity  of  the  corners,  neglected  in 
this  first  approach,  is  currently  being  introduced  into  the  solution.  The  pin 
joint  problem  has  been  carried  out  by  Oplinger  and  Gandhi [6],  and  initial  re¬ 
sults  are  shown  in  Figure  17.  The  solution  here  is  for  a  "rigid"  pin  in  an 
infinite,  homogeneous  but  anisotropic  plate  with  perfect  fit.  A  computer  pro¬ 
gram  has  been  completed.  The  stress  levels  are  given  for  an  applied  pin  load 
of  100  lbs.  per  inch.  The  stress  distribution  is  given  for  the  hole  boundary, 
i.e.,  at  y  =  a,  and  also  for  a  circular  path  at  y  =  1.5a.  To  account  for  the 
presence  of  the  laminae,  a  separate  analysis  is  underway  as  implied  above 
treating  the  physical  problem  shown  in  the  upper  left  of  Figure  15.  Here  an 
infinite  layered  plate  with  a  hole  subject  to  tension  is  being  analyzed.  The 
thickness  variations  associated  with  the  change  in  lamina  modulus  with  each 
layer  are  being  treated.  Results  of  this  solution  will  then  be  integrated 
with  the  rigid  pin  solution  to  obtain  the  total  solution.  The  results  will 
also  be  extended  to  account  for  finite  width  specimens. 

One  of  the  principal  areas  of  potentially  fruitful  exploitation  of  frac¬ 
ture  mechanics  is  in  the  prediction  of  residual  strength  and  crack  arrest 
properties  of  critical  structural  components  of  military  vehicles  (aircraft, 
etc.)  where  those  components  have  been  damaged  in  a  combat  situation.  Figure 
18  suggests  schematically  one  technique  for  arrest.  We  are  initially  under¬ 
taking  a  study  of  the  structural  mechanics  effects  from  damage  due  to  pene¬ 
tration  by  fragments  and/or  projectiles.  Structural  damage  due  to  blast 
loadings  is  considered  for  study  at  a  later  date. 

Rich  [7]  of  our  staff,  after  surveying  the  current  state-of-the-art  of 
battle  damage  analysis,  has  defined  three  approaches  to  projectile  damage: 
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1)  Empirical  approach  based  upon  residual  strength  tests  of  projectile  damaged 
components;  2)  Statistical  approach  based  upon  reliability  studies  of  residual 
strength  tests;  and  3)  Analytical  approach  based  upon  stress  computations  of 
representative  damage  configurations  plus  appropriate  failure  criteria.  In 
the  first  two  approaches  a  measure  of  damage  is  required  to  correlate  results; 
one  such  possibility  is  shown  as  an  "effective"  damage  crack  length  in  Figure 
19.  To  evaluate  all  three  approaches,  we  are  conducting  residual  strength 
tests  on:  a)  controlled,  simulated  battle  damage  components;  and  b)  actual 
ballistically  damaged  components  as  suggested  in  Figure  20.  The  eventual 
product  of  these  studies  will  be  residual  strength  curves  for  the  best  measure¬ 
ment  of  battle  damage.  These  curves  should  aid  both  the  designer  of  future 
equipment  in  selecting  designs  of  reduced  vulnerability  and  the  field  com¬ 
mander  in  assessing  the  operational  limitations  on  current  equipment. 

With  the  exception  of  the  immediately  preceding  discussion  of  residual 
strength,  I  have  been  focusing  attention  on  heterogeneous  materials  and  their 
fracture  behavior.  Since  analysis,  taking  this  heterogeneity  into  account, 
in  any  depth  is  extremely  difficult  and  cumbersome,  considerable  effort  has 
been  devoted  by  the  staff  here  at  AMMRC  to  an  alternative  of  considering  the 
problems  of  homogeneous  but  anisotropic  materials.  The  test  of  the  adequacy 
of  such  an  approach  would  eventually  rest  with  experimental  substantiation. 

Bowie [8]  has  pioneered  in  the  application  of  complex  variable  and  con¬ 
formal  mapping  to  crack  problems.  Basically,  this  concept  involves  the  in¬ 
verse  mapping  of  the  crack  and  its  environment  into  the  unit  circle  and  its 
environment  and  the  finding  of  a  suitable  stress  function  usually  expressed 
as  a  single  Laurent  series.  He  has  recently  modified  this  approach,  devel¬ 
oping  a  technique  which  retains  the  character  of  the  exact  solution  near  the 
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crack  tip  but  which  incorporates  the  numerical  simplifications  associated 
with  the  collocation  technique.  In  this  Modified  Mapping  Collocation  (MMC) 
technique  he  uses  a  relatively  simple  mapping  function  to  detail  only  the 
crack  tip  region,  but  then  enforce*,  satisfaction  of  boundary  conditions  in  a 
collocation  sense.  Figure  21  shows  a  selection  of  configurations  which  Bowie 
has  treated  using  this  technique.  A  more  detailed  representation  of  results 
for  the  specific  case  of  non-uniform  periodic  edge  cracks  is  shown  in  Figure 
22.  This  shows  the  stress  intensity,  normalized,  as  a  function  of  crack 
length  ratios  with  crack  "spacing"  as  a  parameter.  This  problem  was  initial¬ 
ly  motivated  by  the  gun  tube  problem  wherein,  due  to  the  periodic  bore  rifling 
and  the  associated  stress  concentration  factor,  a  periodic  array  of  radial 
cracks  generally  form.  It  is  important  to  know  the  pertinent  stress  intensity 
factors  and  particularly  the  interaction  effects.  The  configuration  of  Figure 
22  was  an  interim  step  toward  this  solution  as  a  development  of  the  cylin¬ 
drical  configuration.  Concurrent  with  this  exercise,  however,  Bowie  and 
Freese  have  in  fact  carried  out  the  case  of  the  finite  wall  "ring"  with  radial 
cracks [9] . 

One  of  the  major  contributions  associated  with  the  MMC  method  is  that  it 
is  appropriate  not  only  for  isotropic  materials  but  also  for  anisotropic 
materials.  It  is  a  highly  flexible  technique.  In  fact,  it  has  been  pro¬ 
grammed  computationally  ("ORTHO")  so  that  wide  classes  of  geometric  config¬ 
urations  are  readily  accessible.  Figure  23  shows  the  types  of  geometries 
which  are  treated  by  "ORTHO".  The  configurations  on  the  left  and  right  col¬ 
umns  are  intended  to  be  rectangular;  they  are  not  restricted  to  square  panels. 
Note  the  eccentricity  of  the  cutouts  or  cracks  and  the  various  other  plate 
configurations,  i.e.,  circular  and/or  elliptical.  And  of  course  these  can 
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be  treated,  as  stated,  either  as  isotropic  or  anisotropic  plates. 

It  is  noted  also  that  the  computational  program  "ORTHO"  has  been  developed 
with  a  high  emphasis  on  user  orientation.  "ORTHO"  has  been  prepared  by  Freese 
so  that  it  can  be  used  by  the  engineer  who  has  no  -  repeat  NO  -  computer  or 
analytical  expertise  to  get  stress  intensity  factors  or  stress  distributions 
in  these  wide  classes  of  problems.  Typical  training  of  a  "user"  of  ORTHO  en¬ 
tails  no  more  than  5-10  minutes  of  instruction  and  the  punching  of  two 
IBM  cards  if  the  material  is  isotropic  and  at  most  three  IBM  cards  for 
anisotropic  materials . 

There  arc,  however,  limitations  to  this  MMC  method,  and  one  such  limita¬ 
tion  is  best  discussed  in  terms  of  the  example  of  Figure  24.  Here  are 
plotted  the  stress  intensity  vs  crack  length  for  a  rectangular  panel  with 
edge  cracks:  a)  with  free  ends  under  uniform  tension;  and  b)  with  constrainet 
ends  under  tension.  Note  that  for  the  free  end  case,  the  stress  intensity 
factor  H  (normalized)  is  relatively  insensitive  to  the  aspect  ratio  w/h  once 
w/h  >_  3.  Thus  the  greatest  ratio  one  would  have  to  consider  in  this  case 
is  then  w/h  =  3.  Such  a  relatively  short  panel  i  .-ily  tractable  using  a 
single  power  series  expansion  to  represent  the  stress  function  in  the  entire 
region  of  the  panel.  On  the  other  hand,  for  the  constrained  end  case,  note 
the  sensitivity  of  the  solution  even  at  aspect  ratios  as  high  as  w/h  =  16. 
Though  intuitively  one  would  expect  the  solution  for  large  w/h's  to  approach 
that  of  the  infinite  plates  with  free  ends  represented  by  the  top  curve 
w/h  =  3,  the  solution  has  in  fact  not  yet  achieved  that  level  even  at  w/h  = 

16.  Solutions  for  large  w/h  are  not  very  tractable  using  a  single  power 
series  expansion  to  represent  the  stress  function.  Bowie  et  al  [10]  have 
therefore  attacked  this  difficulty  by  considering  such  large  or  awkward 


configurations  as  "partitioned"  into  smaller  panels,  and  he  then  associates 
with  each  of  these  partitions  a  separate  power  expansion  to  represent  the 
stress  function  in  the  respective  areas.  The  solutions  in  each  partition 
are  then  "stitched"  together  at  the  interfaces  under  imposed  compatibility 
requirements  in  a  collocation  sense.  Figure  25  shows  the  specific  parti¬ 
tioning  plan  used  in  the  present  constrained  ends  problem.  The  symmetry 
of  the  problem  was  enforced  by  control  of  the  coefficients  and  exponents  of 
the  Laurent  series  in  I.  Expansions  of  power  series  were  also  taken  about 
the  centers  of  regions  of  II  and  III.  Compatibility  of  the  solutions  at  the 
dashed  lines  was  enforced.  Using  this  approach,  the  solutions  plotted  in 
Figure  24  were  readily  obtained.  Currently,  Bowie  is  exploring  the  question 
of  tradeoff  between  using  a  greater  number  of  partitions  vs  phasing  in  of 
finite  elements  in  order  to  optimize  the  size  of  the  numerical  system. 

Bowie  and  Freese[ll]  and  Lakshmikantham[12]  were  among  the  first  investigators 
to  recognize  the  marked  effect  of  orthotropy  in  the  crack  problem.  Figure  26 
from  Bowie  and  Freese [11]  portrays  some  results  of  stress  intensity  factors 
vs  crack  length  in  an  orthotropic  finite  plate  over  a  range  of  aspect  ratios 
c/b.  Figure  27  taken  from  Lakshmikantham[12]  more  dramatically  shows  the 
effect  of  orthotropy  on  stress  intensity  factors.  It  is  of  particular  in¬ 
terest  to  note  the  considerable  difference  between  the  orthotropic  and  the 
isotropic  solution.  Only  at  small  crack  lengths  (A)  or  equivalently  for 
infinite  sheets  do  the  results  merge.  This  perhaps  accounts  for  the  earlier 
misconception  that  material  orthotropy  did  not  influence  stress  intensity 
levels . 

Gandhi [13]  treated  the  more  general  case  of  arbitrary  orientations  of 
the  crack  and  material  orthotropy  relative  to  the  loading  direction.  Figure 
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28  shows  some  representative  results  for  the  special  case  of  a  square  plate 
with  a  crack  at  45°  and  portrays  the  influence  of  material  orientation  on 
the  normalized  stress  intensity  factor  for  both  the  Mode  I  (crack  opening) 
and  Mode  II  (in-plane  shear)  conditions. 

Lakshmikantham[12]  has  examined  also  the  effect  of  stiffness  in  an  in¬ 
finite  strip  of  orthotropic  material  with  cracks.  Results  given  in  Figure 

29  show  the  behavior  with  relatively  low  axial  stiffness,  intermediate  and 
high  axial  stiffness.  It  is  significant  to  note  that  for  the  low  stiffness 
case  the  stress  intensity  curves  for  both  isotropic  and  orthotropic  panels 
display  an  ever  increasing  slope  suggesting  a  rather  unstable  situation;  i.e., 
when  a  crack  starts  to  grow,  it  probably  continues  to  expand  in  a  catastrophic 
fashion.  On  the  other  hand,  as  the  stiffness  is  increased,  the  slopes  tend  to 
level  off,  and  in  fact  for  the  higher  stiffnesses  shown  the  curves  take  on 
negative  slopes  suggesting  a  more  stable  behavior,  and  in  fact  suggest  the 
feasibility  of  crack  arrest.  In  the  case  of  the  highest  stiffness  the  ortho¬ 
tropic  panel  exaggerates  this  negative  slope  even  more  than  the  isotropic 
case  suggesting  its  more  favorable  potential  in  a  crack  environment.  Analysis 
of  wide  plates  leads  to  the  same  general  results[14]. 

Thus  far  in  the  discussion.  I’ve  presented  only  our  static  considera¬ 
tions.  We  are  also  embarking  on  a  serious  effort  to  examine  dynamic  aspects 
of  fractures,  particularly  in  the  low  microsecond  regime.  Jones  has 
recently  completed  the  installation  and  preliminary  testing  of  a  4-inch 
diameter  light  gas  gun  for  use  in  plate  slap  studies.  Figure  30  shows  the 
facility  in  place.  The  breech  assembly  is  to  the  left,  the  launch  tube  is 
approximately  26  feet  long,  feeds  into  t'  .arget  area,  and  thence  into  the 
large  target  recovery  chamber  clearly  visible  at  the  right  end.  Figure  31 
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is  a  schematic  view  of  the  system  (oriented  in  juxtaposition  to  the  previous 
figure) .  This  system  is  capable  of  achieving  impact  pressures  up  to  570 
kilobars  and  impact  velocities  up  to  4000  feet  per  second  dependent  upon  the 
specimen  mass.  The  curves  at  the  upper  half  of  the  figure  suggest  the  type 
of  data  we  expect  to  obtain;  specifically,  equations  of  state  data  including 
the  Hugoniot  yield  value  and  fracture  stress  in  a  shock  environment,  i.e., 
spallation  will  be  determined.  "Projectiles"  fired  are  in  the  form  of  thin 
discs  supported  during  launch  by  a  sabot  arrangement.  Targets  are  also  disc¬ 
like.  The  curve  at  the  upper  left,  though  identified  as  impact  velocity  vs 
projectile  thickness,  can  be  interpreted  as  impact  stress  and  pulse  duration 
respectively,  i.e.,  by  varying  the  impact  velocity,  the  stress  level  can  be 
controlled,  and  by  varying  the  projectile  thickness,  the  pulse  duration  can  be 
conti'olled. 

Requirements  for  dynamic  fracture  criteria  are  motivated  by  our  need  to 
incorporate  such  data  into  penetration  mechanics  studies.  Computer  codes, 
particularly  HEMP*,  are  being  exploited  by  Mescal 1  in  treating  the  armor  - 
projectile  interaction  problem.  A  typical  computer-simulated  result  is  shown 
in  Figure  32.  This  shows  the  state  of  penetration  some  0.78  micro  seconds 
after  impact.  Since  1'otationa]  symmetry  is  assumed,  the  view  shown  is  a 
radial  plane  about  the  centerline  of  impact.  The  projectile  impacting  from 
tiie  left  at  a  velocity  of  7,000  feet  per  second  was  initially  a  regular 
cylinder  of  steel;  appreciable  deformation  is  clearly  visible.  The  armor 
system  consisted  of  a  front  layer  of  boron  carbide  B^C,  a  relatively  brittle 
ceramic,  and  a  backup  plate  of  aluminum.  Note  the  two  failure  areas  of  the 
forward  ceramic  plate.  In  particular,  the  failure  by  shear  in  this  brittle 

*HEMP  -  developed  by  M.  L.  Wilkins,  Lawrence  Livermcre  Labs. 
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ceramic  was  predicated  on  a  maximum  plastic  effective  strain.  This  may  ap¬ 
pear  puzzling,  i.e.,  that  plastic  -  strain  is  tolerated  in  such  a  material. 

On  the  other  hand,  it  must  be  kept  in  mind  that  this  is  possible  only  be¬ 
cause  of  the  excessively  high  hydrostatic  pressure  existing  in  the  region 
of  impact.  Under  such  ambient  pressures  brittle  materials  may  deform 
plastically.  The  actual  strain  to  fracture  was  in  fact  taken  to  be  a  func¬ 
tion  of  the  pressure/effective  stress  ratio.  The  exact  nature  of  this  func¬ 
tion  is  not  known  in  detail;  rather  only  a  few  points  are  well  characterized. 
More  detailed  information  is  needed.  We  will  be  seeking  some  data  to  fill 
this  gap. 

In  the  second  area  of  failure,  identified  as  tension  failure,  it  was 
determined  that  this  region  was  subjected  to  tensile  stresses  in  excess  of 
30  kilobars.  This  is  vastly  in  excess  of  the  anticipated  fracture  strength 
of  the  material,  but  no  hard  data  are  available  on  these  values  -  so  once 
again  it  is  evident  that  fracture  data  -  characterization  of  these  rela¬ 
tively  brittle  materials  under  short  duration-high  intensity  tension  pulses  - 
is  required  if  we  are  to  reliably  predict  performance  in  the  penetration 
environment . 

III.  FATIGUE 


Our  studies  in  fatigue  have  been  somewhat  more  limited  than  the  frac¬ 
ture  investigations.  We  have  one  team  of  investigators  exploring  the  crack 
initiation  and  another  the  crack  propagation  stages  of  fatigue.  The  emphasis 
currently  in  both  these  efforts  relates  to  behavior  of  metals.  Composite 
behavior  under  fatigue  is  only  beginning  to  receive  more  of  our  attention. 

The  crack  initiation  study  led  by  Papimo  is  aimed  at  acquiring  some 


13 


quantitative  observations  of  crack  initiation  mechanisms  and  then  at  the 
establishment  of  pragmatic  criteria  for  crack  initiation  which  can  be  inte¬ 
grated  with  cyclic  crack  propagation  and  ultimate  fracture  criteria  to  pre¬ 
dict  failure  of  components. 

As  a  preliminary  step  aimed  at  defining  and  quantifying  "initiation", 
Papimo  et  al  [15]  [16]  have  developed  an  automatic  system  to  take  flash 
photomicrographs  at  roots  of  notches  in  flat,  doubly-notched  specimens  during 
low  cycle  fatigue  tests.  A  schematic  view  of  the  notch  region  of  the  specimen 
is  shown  in  Figure  33.  A  wide  range  of  theoretical  elastic  stress  concentra¬ 
tion  factors  is  readily  obtained  with  this  configuration.  The  present 
study  has  been  limited  to  date  to  a  range  of  K  from  1.5  to  5.0.  Figure  34 
shows  the  arrangement  of  apparatus  for  the  photomicrography.  It  is  noted,  of 
course,  that  with  the  doubly  notched  specimen  of  Figure  33  and  the  dual 
photographic  assembly,  each  test  constitutes  two  sets  of  data,  one  for  each 
notch. 

The  photographic  apparatus  is  synchronized  with  the  loading  machine  to 
take  photomicrographs  on  35mm  film  at  preprogrammed  cyclic  intervals  without 
arresting  the  cyclic  loading  procedure.  This  is  continued  until  a  well  de¬ 
fined  crack  has  formed.  The  developed  film  is  then  examined  and  the  crack  is 
traced  backward  in  time  frame  by  frame.  A  crack  was  judged  to  have  just 
initiated  when  it  was  visible  in  one  frame  but  not  in  the  frame  just  pre¬ 
ceding.  Resolution  of  cracks  with  a  surface  length  of  .001  inch  and  an  open¬ 
ing  of  .0001  inch  was  easily  effected  using  projected  images.  Repetition 
rates  for  recording  photos  was  generally  100  cycles,  though  a  limited  number 
of  tests  at  veiy  low  cycle  life  were  photographed  at  50  cycle  intervals. 
Furthermore,  the  photomicrographs  were  taken  at  maximum  load  in  the  cycle,  thus 
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making  the  cracks  more  visible.  These  techniques  permitted  confident  identi¬ 
fication  of  crack  initiation  to  within  ±50  or  ±25  cycles. 

Figure  35  shows  a  selected  sample  of  photomicrographs  of  the  notch  base 
after  various  number  of  cycles.  These  were  taken  from  a  long  sequence  of 
films.  The  particular  specimen  involved  had  an  initial  elastic  stress  con¬ 
centration  factor  of  of  4,  and  the  figure  illustrates  one  form  of  crack¬ 
ing  behavior  where  a  number  of  individual  cracks  initially  form  and  then 
coalesced  into  one  dominant  crack  and  a  few  dormant  cracks.  The  initiating 
cracks  are  more  clearly  visible  in  the  negatives.  The  circle  outlines  the 
site  of  the  first  crack.  Generally,  though  not  consistently,  the  initial 
cracks  developed  near  the  mid  thickness.  Furthermore,  the  growth  pattern 
was  not  always  one  of  several  initiation  cracks  coalescing  into  one  dominant 
crack.  More  generally  -  and  this  depends  not  only  on  the  materials  but  also 
on  the  initial  stress  concentration  factor,  and  the  load  (and  quite  probably 
on  the  gradient  of  stress)  -  only  some  of  the  initial  cracks  coalesced  to 
form  the  dominant  crack;  others  remained  isolated  and  dormant.  Figure  36 
shows  a  typical  count  of  the  cracks  visible  at  the  notch  base  as  a  function  of 
cycles.  It  is  noted  that  both  notches  developed  first  cracks  at  about 
1600-1700  cycles;  subsequently,  as  many  as  six  independent  cracks  were  ob¬ 
served  on  one  side  (four  on  the  other),  eventually  a  dominant  crack  became 
evident  (at  about  5000  cycles)  though  the  dominant  crack  was  not  the  only 
crack  present.  The  indicated  points  correspond  to  the  cycle  count  at  which 
the  crack  intercepted  the  face  of  the  specimens .  In  this  particular  case  one 
crack  initiated  at  the  edge;  the  other  required  some  3000  cycles  before  reaching 
the  edge.  This,  of  course,  highlights  the  difficulty  associated  with  interpre¬ 
tations  based  solely  on  observations  made  on  the  faces  of  notched  specimens. 
Initial  testing  was  carried  out  with  zero  -  tension  loading.  Additional  tests 


are  being  done  over  the  range  of  compression  to  tension. 

The  limited  tests  carried  out  to  date  have  been  on  an  AISI  4340  steel. 
Selection  was  not  based  on  an  inherent  interest  in  this  particular  material 
but  rather  it  was  desired  to  use,  as  a  vehicle  for  study  of  and  development 
of  pragmatic  design  criteria,  a  well  characterized  and  structurally  signifi¬ 
cant  material.  Other  more  advanced  materials  such  as  Transformation  Induced 
Plasticity  (TRIP)  steel,  for  example,  are,  however,  to  be  explored. 

As  an  adjunct  to  these  studies  and  in  anticipation  of  the  more  detailed 
continuum  treatment  of  the  elastic-plastic  behavior  of  the  notched  material, 
Papirno  and  Campo[17]  have  developed  a  technique  for  extending  the  strain  range 
of  conventional  strain  gages.  Figure  37  shows  the  arrangement  used.  With  this 
procedure,  cyclic  strains  of  the  order  1-3%  are  conveniently  measured  and 
monitored.  Satisfactory  results  have  been  achieved  for  as  many  as  5000  cycles 
at  strain  levels  „  1.5%. 

In  exploring  crack  propagation  behavior,  Matthews  et  al [18]  have  re¬ 
cently  reported  effects  of  sudden  or  step  changes  in  the  stress  intensity 
factor  on  propagation  rates.  They  tested  sheet  specimens  of  AISI  4340  steel 
quenched  and  tempered  to  an  ultimate  tensile  strength  of  1S1  ksi.  The  speci¬ 
men  configuration  along  with  some  results  are  sketched  in  Figure  38.  Speci¬ 
mens  0.157  and  0.025  inches  in  thickness  were  cycled  (tension  -  tension; 
K]iin/*Snax  =  0.1)  in  a  dry  air  atmosphere  using  a  closed  loop  hydraulic  testing 
machine.  The  load  was  reduced  incrementally  at  crack  extension  intervals  of 
0.025  inches  to  maintain  a  quasi-constant  AK  loading  to  within  ±1.3%  of  the 
desired  level.  This  constant  AK  was  maintained  until  the  growth  rate  reached 
a  steady  state.  At  this  point  the  load  level  was  raised  to  obtain  a  new 
higher  AK  level.  The  cyclic  loading  was  reapplied  and  the  crack  growth  rate 
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again  observed.  During  this  second  stage,  also,  the  load  was  adjusted  peri¬ 
odically  to  maintain  the  constancy  of  the  augmented  AK  level  as  the  crack 
progressed.  The  crack  length  was  monitored  at  frequent  cyclic  intervals  using 
a  traveling  microscope  on  both  sides  of  the  specimen.  Typical  results  shown 
in  Figure  36  show  a  rather  significant  transient  augmentation  of  the  crack 
growth  rate  well  beyond  that  steady  state  level  which  could  be  expected  based 
solely  upon  the  second  stage  stress  intensity  factor.  This  is  an  unconserva¬ 
tive  trend  and  could  have  serious  consequences  in  critical  designs  if  not 
recognized  and  taken  into  consideration.  We  are  attempting  to  formulate 
models  which  will  adequately  and  conveniently  account  for  this  behavior. 
Studies  will  also  include  observation  of  behavior  when  the  second  stage  AK 
range  is  decreased  from  the  level  of  the  first  stage. 

In  a  separate  but  related  effort  to  try  to  provide  a  useful  tool  for 
prediction  of  steady  state  crack  propagation  rates  over  the  wide  spectrum  of 
AK  values  from  threshhold  to  critical  levels,  Bluhm  et  al[19]  have  formulated 
the  following  tentative  and  empirical  relation  which  does  in  fact  account  for, 
at  least  qualitatively,  transitions  in  crack  propagation  behavior  over  these 
extremes  of  AK  levels. 
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where  C,  m,  n  are  constants, 

AKq  is  the  threshhold  level  of  the  stress  intensity  factor  range 
corresponding  to  Kmax, 

K  is  the  maximum  cyclic  value  of  the  stress  intensity  factor, 
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Kc*  is  an  appropriate  fracture  toughness  value  dependent  upon  , 

plate  thickness,  yield  strength,  etc.  and  based  upon  the  model  for 
thickness  effect  described  by  Bluhm[20]. 

We  have  not,  however,  had  the  opportunity  yet  to  experimentally  verify  the 
range  of  validity  of  this  relation.  The  form  of  this  relation  is  of  partic¬ 
ular  interest  in  that  it  displays  a  non-uniqueness  of  da/dN  for  a  given  AK; 
because  the  thickness  effect  is  explicitly  reflected  in  Kc*,  the  da/dN  curve 
takes  on  the  characteristic  behavior  sketched  in  Figure  39. 

Beeuwkes[21]  has  examined  the  more  restricted  crack  propagation  "law" 
defined  in  the  general  form 


§  =  <^2 


(3) 


— 2  — 2 

Replacing  AK  by  its  stress  equivalent  Ao  a  (for  a  small  crack),  this  equa 
tion  can  be  written  as  follows* 


da 

a 


=  C  Ao  dN 


and  by  integration, 


In  a  =  In  a  +  C  Ao  (N  - 


(4) 


(5) 


*Finite  width  corrections  using  the  Feddersen  secant  correction  factor  were 
actually  considered  by  Beeuwke|  and  led  to  Equation  fa)  below  instead  of 
Equation  (4)  above  where  is  the  cosine  integral.  (Nevertheless  one 

is  led  to  the  same  conclusions.) 

ci(r}  =  C0  +  Ci^2  (N  -  V  (a) 
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This  form  of  a  versus  N  leads  to  a  straight  line  prediction  for  In  a  versus 
N.  In  the  plotting  of  extensive  data  Beeuwkes  noted,  however,  that  the  data 
in  fact  were  better  represented  by  a  series  of  intersecting  straight  lines 
with  each  successive  segment  having  approximately  twice  the  slope  of  the  pre¬ 
ceding  segment.  Typical  data  are  shown  in  Figure  40.  These  transitions 
were  not  always  associated  with  observable  gross  changes  in  fracture  surface 
appearance.  This,  however,  is  being  explored  further. 

When  da/dN  is  appropriately  determined  graphically  from  this  data  one 
is  consistently  led  to  a  discontinuous  log-log  plot  of  da/dN  versus  AK  as 
shown  in  Figure  41.  Of  particular  interest  is  Beeuwkes*  observation  that 
the  slopes  of  these  disjointed  segments  are  statistically  identical  and  ap¬ 
proximately  2  which  is  consistent  with  the  exponent  of  Equation  3.  However, 
it  is  further  noted  that  the  mean  "curve"  going  through  all  the  segments  has 
the  slope  4  which  has  frequently  been  identified  as  the  appropriate  power. 
Some  100-200  crack  propagation  curves  from  a  dozen  or  so  investigators  have 
been  analyzed  by  this  treatment  and  almost  without  exception  behave  in  iden¬ 
tical  fashion. 

IV.  PERIPHERAL  STUDIES 

Contributing  to  and  complementing  both  the  fracture  and  fatigue  studies, 
is  a  program  being  carried  out  by  Lenoe  aimed  at  defining  constitutive 
relations  for  various  classes  of  structurally  important  materials.  He  has 
initiated  a  series  of  creep  and  relaxation  tests  on  unidirectional  fiber  re¬ 
inforced  composites.  To  date  the  materials  studied  have  been  S  glass,  boron, 
Thornel  50  and  Mod  Mor  II  reinforced  epoxies  over  the  temperature  range  from 
room  temperature  to  400°F.  Representative  creep  and  relaxation  data  are 


shown  in  Figures  40  and  41.  Additional  materials  to  be  investigated  will 
include  E  glass  in  various  epoxy  resins,  PRD-49,  as  well  as  metal  and  ceramic 
matrix  composites.  Studies  will  also  be  extended  to  include  angle  ply  con¬ 
figurations.  Testing  will  be  in  longitudinal  and  transverse  tension,  torsion 
and  in-plane  shear. 

A  major  objective  of  this  investigation,  aside  from  the  data  collection 
aspect,  is  to  examine  various  theoretical  approaches  to  creep  prediction 
from  relaxation  response  (or  vice  versa)  and  to  provide  techniques  for  esti¬ 
mating  angle  ply  behavior  from  known  behavior  of  unidirectionally  reinforced 
composites.  The  most  appropriate  of  these  approaches  might  then  be  recom¬ 
mended  as  an  engineering  analysis  procedure. 

Additionally,  Chou  has  been  examining  the  dynamic  behavior  of  several 

polymers;  specifically  the  stress-strain  curves  at  strain  rates  up  to  approxi- 
3 

mately  10  per  second,  particularly  in  compression,  were  determined.  Prelim¬ 
inary  results  revealed  that  one  polymer,  polypropylene,  displayed  rather  sur¬ 
prisingly  a  maximum  load  feature;  this  is  usually  observed  only  in  tension 
specimens.  In  attempting  to  account  for  this  behavior  and  since  this  material 
was  known  not  to  have  a  phase  transition  at  the  nominal  stresses  and  temper¬ 
ature  of  the  test,  Chou,  anticipating  a  possibly  significant  temperature  rise, 
instrumented  a  number  of  specimens  with  imbedded  thermocouples  and  repeated 
the  tests. 

Results  of  several  representative  compression  tests  on  the  polypropylene 
are  shown  in  Figure  44.  Three  stress- strain  curves  are  shown  for  strain 

_3 

rates  from  2  x  10  /sec  to  1530/sec.  Only  the  curve  for  the  high  strain 
rate  shows  this  maximum  load  characteristic.  In  the  same  figure  are  plotted 
the  recorded  temperatures.  It  is  noted  that  at  the  two  higher  strain  rates, 
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a  marked  knee  in  the  temperature,  stress  curve  is  evident.  At  the  highest 
strain  rate  (e  =  1530/sec)  this  knee  occjrs  at  the  maximum  load  point.  We 
are  not  sure  yet  of  the  mechanisms  interacting  to  generate  the  temperature 
rise  and/or  the  peaking  behavior  of  the  stress-strain  curve. 

Our  interest  here  is  motivated  by  the  practical  implication  with  respect 
to  the  consequences  of  such  behavior  in  these  materials  if  used  in  personnel 
armor  or  more  generally  as  a  matrix  material  in  composites. 

V.  CONCLUSION 


Summing  then,  I  have  briefly  described  a  sampling  of  the  experimental 
and  analytical  studies  at  AMMRC  relating  to  fracture  and  fatigue  and  have 
even  more  briefly  made  reference  to  some  of  our  activities  in  constitutive 
relationships.  Time  limitations  have  precluded  more  detailed  description 
of  the  foregoing  or  reference  to  other  relevant  studies,  such  as  the  follow 
ing: 


Beeuwkes 

Baratta 

Lenoe 


on  Fracture  Criteria  and  Crack  Growth  Law. 
on  Design  of  Specimens  for  Fracture  Testing  of 
Brittle  Materials. 

on  Integration  of  Specimen  Configurations  and 

Fracture  Considerations  in  Reliability  Considera- 

'V  • 

tions . 


Papirno  on 

Matthews  on 

Interested  personnel 
further  elaboration. 


Elastic-Plastic  Notch  Behavior. 

his  "Collection  of  Plane  Strain  Fracture  Toughness 

Data". 

are  invited  to  contact  the  investigators  directly  for 
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CONSTITUTIVE  RELATIONS 

FATIGUE 

FRACTURE 

WEAR  AND  EROSION 
STRUCTURAL  DYNAMICS 
STRUCTURAL  INSTABILITY 
STRESS  AND  DEFLECTION 
STRUCTURAL  OPTIMIZATION 

Figure  1.  Solid  Mechanics  Responsibility  at  AMMRC 


25 


Compliance 


D> 


s/M;-. 


1  9-066-801  /AMC-72 


Figure  3.  Compact  Tension  Specimens  for  Composite  Fracture  Studies 
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Figure  4.  Compliance  Data  for  Unidirectional  S-Glass/Epoxy 
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Fiber  Direction 
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a.  Unidirectional  fiber  specimen. 


FIBER  ORIENTATION 
(at  surface) 


Figure  7.  Effect  of  Fiber  Orientation  on  Crack  Pattern 
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Figure  9.  Deformation  and  Interlarniner  Shear  Distributions  in  Angle  Ply  Composites 
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Figure  10.  Moire  Patterns  in  ±20  Degree  Boron  Epoxy  Tensile  Specimen 
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Figure  12.  Distribution  of  Interlaminar  Shear  Stress  t%z  at  Free  Edge  (y=b)  of  Angle  Ply  Strip  under  Tension. 
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YARN  REINFORCEMENT 
(graphite,  glass,  etc.) 
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Figure  13.  Effect  of  Filament  Structure  on  Deformation  and  Interlaminar  Shear  Stress  at  Edge  y=b 
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STACKING  SEQUENCES 


C  (CLUSTERED) 


4 1+9 
SI -6 
4 1+0 


16 

Layers 


A  (ALTERNATING) 

::|  * 

Layers 


PLY 

ANGLE 


BORON 

EPOXY 


GRAPHITE 

EPOXY 


SEQUENCE  SEQUENCE 

A  C _ A  C 

XX  XX 


X  X 
X  X 


Figure  14.  Matrix  of  Experiments 
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Loads 


Figure  15.  Typical  Layered  Structural  Configuration  under  Study  in  Connection  with  Joints 


Figure  18.  Laminated  Stri 
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LD:  MAXIMUM  LATERAL  DAMAGE  DIMENSION. 

TLD:  LATERAL  DAMAGE  DIMENSION  TRANSVERSE 
TO  MAIN  LOADING  DIRECTION. 


Figure  19.  Illustration  Showing  Effective  Damage  "Length" 
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Figure  21 .  Representative  Configurations  Treated  by  Modified  Mapping  •  Collocation  Method 
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Figure  24.  Stress  Intensity  in  Rectangular  Panel  with  Edge  Cracks  (Typical  Results) 


Figure  26.  Stress  Intensity 
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Figure  28.  Variation  of  Normalized  Stress  Intensity  Factor  with  Material  Orientation  Angle  0, 
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Figure  30.  Photograph  of  4-inch  Plate  Accelerator  Facility 
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Schematic  of  4-inch  Plate  Accelerator  Facility 
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Figure  32.  Computer  Simulation  of  Penetration  in  Composite  Armor  System 
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Figure  35.  Crack  Initiation  Progression 
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Figure  37.  Strain  Gage/Pad  Configuration 
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Figure  38.  Transient  Crack  Growth  Behavior  after  Step  Change  in  Loading 
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Figure  40.  Typically  Observed  Crack  Length  versus  Number  of  Cycles 
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